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Abstract 

It has been known for some time that the well established LMA solution to 
the observed solar neutrino deficit fails to predict a flat energy spectrum for Su- 
perKamiokande as opposed to what the data indicates. It also leads to a Chlorine rate 
which appears to be too high as compared to the data. We investigate the possible 
solution to these inconsistencies with non standard neutrino interactions, assuming 
that they come as extra contributions to the v a vp and v a e vertices that affect both 
the propagation of neutrinos through solar matter and their detection. We find that, 
among the many possibilities for non standard couplings, only one of them leads to 
a flat SuperKamiokande spectral rate in better agreement with the data and predicts 
a Chlorine rate within la of the observed one, while keeping all other predictions 
accurate. 
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1 Introduction 



Neutrino non-standard interactions (NSI) have been introduced long ago [H [2] to account 
for a possible alternative solution to the solar neutrino problem. Since then a great deal 
of effort has been dedicated to study its possible consequences. To this end possible NSI 
signatures in neutrino processes have been investigated, models for neutrino NSI have been 
developed and bounds have been derived [3l SI El d [HI [91 [101 [TU [121 UHl HS1 [16l HZ] - 
Specific investigations of NSI in matter have also been performed within the context of 
supernova [TS] and solar neutrinos [151 1201 12T1 1221 1231 123] . 

Although LMA is generally accepted as the dominant solution to the solar neutrino 
problem [251 EE], not only its robustness has been challenged by NSI, as it can shift the 
LMA solution to the dark side region of parameter space [22J , but also some inconsistencies 
remain regarding its agreement with the data [271 EH] . In fact, while the SuperKamiokande 
(SK) energy spectrum appears to be flat [291 [30], th e LMA prediction shows a clear negative 
slope in the same energy range. With the expected improvement in the trigger efficiency 
for threshold electron energies as low as 3 MeV to be reached in the near future [21], such 
a disagreement, if it persists, may become critical. Moreover the LMA solution predicts an 
event rate for the CI experiment |32j which is 2a above the observed one [33J. These are 
motivations to consider 'beyond LMA' solar neutrino solutions in which NSI may play a 
subdominant, although important role. 

In order for NSI to be detectable and therefore relevant in physical processes, the char- 
acteristic scale of the new physics must not be too much higher than the scale of the physics 
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giving rise to the Standard Model interactions, Aew — G F . Possible realisations are one 
loop radiative models of Majorana neutrino mass [31], supersymmetric SO(10) with broken 
D-parity [35], the inverse seesaw in a supersymmetry context [36] or triplet seesaw models 
[37] . Since the scale at which the new interaction arises is supposed to be not too far from the 
electroweak scale, its coupling may be parameterised by Gpe where e ~ K 2 EW / K 2 NP ~ 10 -2 
for d=6 or e ~ A FW /Aj^ p ~ 10~ 4 for d=8 operators respectively. For type I seesaw, NSI 
are of course negligible. 

In this paper we will be concerned with NSI both at the level of propagation through 
solar matter and at the level of detection. Matter NSI are defined through the addition of 
an effective operator to the Lagrangian density 

C% SI = -2V2G F s^[h"Pf] Ki.PlM (1) 

where / = e,u,d, P denotes the projection operator for left and right chirality and £^ 
parameterizes the deviation from the standard interactions. At present there is no evidence 
at all of such operators generated at a scale Ajvp, hence the variety of theoretical models 
for the physics accessible to the LHC. 

There are several ways to introduce NSI. For instance in fermionic seesaw models, once 
the heavy fermions (singlets or triplets) are integrated out, modified couplings of leptons to 
gauge bosons are obtained in the form of a non-unitary leptonic matrix. The strong bounds 
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on the deviation from the unitarity of this matrix constrain these NSI to be < O(10~ 3 ) [13J. 
Alternatively NSI can be generated by other new physics above the electroweak scale not 
related to neutrino masses. As a consequence an SU(2) gauge invariant formulation of NSI 
is required, since any gauge theory beyond the standard model must necessarily respect its 
gauge symmetry. Strong bounds from four charged fermionic processes [38] and electroweak 
precision tests requiring fine tunings imply that possibilities are limited for such scenarios 
[Hl[39]. Another way to introduce NSI is by assuming that these are extra contributions to 
the vertices v a v$ and u a e. In such a case the parameters e a p describe the deviation from the 
standard model vertices and are treated like the standard interactions. It is possible that 
other effects are present in this case that depend on the nature and number of particles that 
may be introduced in a particular model. We adopt this procedure in the present paper and 
assume these model dependent effects to be negligible. 

The paper is organized as follows: section 2 is devoted to the study of the propagation 
and detection of solar neutrinos. We start by reviewing the derivation of the neutrino 
refraction indices with standard interactions (SI) and their generalization to NSI in order 
to obtain the matter Hamiltonian. The survival and conversion probabilities to and 
v T are then evaluated through the numerical integration of the Schrodinger like equation 
using the Runge-Kutta method and the experimental event rates are obtained. We use the 
reference solar model with high metalicity, BPS08(GS) [40] . In section 3 we investigate the 
influence of the NSI couplings on these rates in order to find whether and how the fits can 
be improved with respect to the LMA ones. We concentrate in particular on the elimination 
of the upturn in the SuperKamiokande spectrum predicted by LMA for energies below 8- 
10 MeV not supported by the data [2Z1 123 1231 EE3| and on the Chlorine rate whose LMA 
prediction exceeds the data by 2a [32, 33j. We find that there are quite limited possibilities 
for the NSI couplings allowing for such fit improvements. Finally in section 4 we draw our 
main conclusions. 

2 Interaction potentials, the Hamiltonian and the rates 

In this section we develop the framework that will be used as the starting point for the anal- 
ysis of the NSI couplings in section 3. To this end we review the derivation of the neutrino 
interaction potentials in solar matter, its generalization to non standard interactions along 
with the corresponding matter Hamiltonian and the event rates. 

2.1 Interaction potentials and the Hamiltonian 

While Ue's propagate through solar matter their interaction with electrons proceeds both 
through charged and neutral currents (CC) and (NC). Recalling that for standard interac- 
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tions (SI) each tree level vertex accounts for a factor 



9l 



cosO 



ir 



(T 3L -2Q f stn 2 6 w ) } (2) 



inserting the W, Z propagators and the electron external lines, one gets for the v e interaction 
potentials 

{V e )cc = G F V2N e , (V e ) NC = G F /V2{-1 + Asin 2 9 w )N e (3) 

where G F is the Fermi constant, G F j\[2 = g1/8 m \v For the interactions with quarks 
only neutral currents are involved and the additivity of the quark-current vertices gives for 
protons 

V p = (V P ) NC = G F /V2(1 - Asin 2 6 w )N e (4) 

and for neutrons 

V n = (V n ) NC = -G F /V2N n . (5) 
Hence the neutrino interaction potential is for standard interactions @ 



V(SI) = V e + V p + V n = GV2N e (l - =V C + V n 



(6) 



with V e = (V e ) cc + (V e ) NC and V c = V e + V p = G F ^2N e . 

In order to introduce NSI we assume that each diagram associated to neutrino propaga- 
tion in matter (i.e. CC and NC currents in v a e~ and NC currents in v a u, v a d scattering) 
is multiplied by a factor s e ^ d F parameterising the deviation from the standard model. So 
we assume that the interaction potential for v a (a = e, //, r) on electrons involves both CC 
and NC giving rise to possible lepton flavour violation: v a for a^e may have CC. So for 
the charged current of v e with electrons we have 

(V e )cc(NSI) = ^r{e%)ccN e (7) 

and for the neutral current 

(V e ) NC (NSI) = + Asin 2 6 w ){e%) NC N e (8) 

where the NSI couplings affecting the CC and NC processes should in principle be distin- 
guished. 

Using equation (2) and the additivity of the quark-current vertices, one gets for the 
neutrino interaction potential with protons 



V P (NSI) - 9 * 



2m 2 w 



u p £ a/3 f^ 0c .uP ^ dp . 

e a/3 2 I 3 °# ~ 3 Q/3 I w 



N e . (9) 



* All expressions are divided by 2 to account for the fact that the medium is unpolarized. 
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Similarly for neutrons 



V n (NSI) 



gi 

2m 2 w 



-uP 

_ pdP I AT 

2 £ a,3 iV «- 



(10) 



In both (9) and (10) only neutral currents are involved. 

Adding (7), (8), (9) and (10) and dividing by 2 one finally gets 



V(NSI) = GV2N e 



{ £ a P p)cC + 



1 2 

-- + -.' 

2 3 



-a/3 



111 



which reduces to eq.(6) in the absence of NSI. 

In the case of the standard interactions, the interaction potentials for v e and u a constitute 
a diagonal matrix because they cannot be responsible for flavour change [eq.(6)]. This may 
occur as a consequence of the vacuum mixing angle (oscillations) [221 EE] or the magnetic 
moment for instance [271 EH]. O n the other hand, in the case of NSI the interaction potentials 
[eq.(ll)] constitute a full matrix in neutrino flavour space. 

In order to obtain the matter Hamiltonian eqs.(6) and (11) must now be added. In the 
flavour basis this is 





( 1 





0\ 




U M = V c \ 











H 




V o 










v ee (NSI) v e ^(NSI) v er (NSI) 
v^ e (NSI) v^(NSI) v^ T (NSI) 
v Te (NSI) v T p(NSI) v TT (NSI) 



(12) 



where in the first term, describing the standard interactions, the additive quantity V n which 
is proportional to the identity, has been removed from the diagonal. In the second term v a p 
(a, (3 — e, fi, t) denote the matrix elements of the interaction potential matrix (11). Finally 
in the mass basis 
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f 









\ 



2E 





2E 



U^HmU 



(13) 



where U is the PMNS matrix [JT] §, E is the neutrino energy and Am?- 



m" with 



i 2 

rrii (i = 1,2,3) the neutrino mass. Upon insertion of this Hamiltonian expression in the 
neutrino evolution equation, the survival (P ee ) and conversion probabilities (P ef i, Per) are 
evaluated using the Runge-Kutta numerical integration. 



§We use the standard parameterization [42] for the U matrix and the central value sinOis = 0.13 claimed 
in ref . [25] . 
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2.2 Neutrino electron scattering detection rates 



For the detection in SuperKamiokande and SNO through v a e — > vp e scattering, the NSI 
information comes in the probabilities and the cross section 



da 2G 2 F m e 
dT = n 



9l + 9r[ 1- 
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where gi,R are the ql,r, couplings modified according to [TT 
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with a = e, /i, r. 

For z/ e both charged and neutral currents are possible, so that 

(9l)u £ = ^ + sin 2 ^v , = sin 2 ^ 

whereas i/^ r only interact through neutral currents, hence 



(15) 



sin 2 9 



w ■ 



(16) 



These expressions are then inserted in the spectral event rate 



R 



■th 

SK,SNO 



dE' e f(E' e} E e ) dE<j>(E) 



da, 



da,, 



da T 



Pee(E)^ + P e »(E)-^ + P eT (E)-^ 



dT 



dT 



dE'J(E' e ,E e ) dE<j>{E) 



da e 
] dT 



(17) 

which will be evaluated in the next section. Here <p(E) denotes the neutrino flux from Boron 
and hep neutrinos, f(E e , E e ) is the energy resolution function for SuperKamiokande and 
SNO [43], HI] and the rest of the notation is standard. 

Notice that whereas the Hamiltonian (13) is symmetric under the interchange 

e a/3 e % for e , U, d 

such is not the case for the detection process [see eq.(14)-(16)]. We finally note that at 
the detection level the NSI couplings e a 'g are considered separately, as clearly seen from 
eqs.(14), whereas at the level of propagation, since the diagrams involved in the interaction 
potentials add up, their sum should instead be considered. 
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3 NSI couplings, probabilities and spectra 



We now perform an investigation of the effect of the NSI couplings e e ^' d = \e e ^ d \e^°^ 
on the neutrino probability and event rates. Our aim is to find those couplings which lead 
to a flat SuperKamiokande spectral rate, thus improving the fit with respect to its LMA 
prediction while keeping the quality of the other solar event rate fits. We first consider 
equal CC and NC couplings, namely (s^)cc = {^%)nc = anc ^ i n a second stage 
{ £ %)cc 7^ {z%)nc- 

For the sake of clarity we will organize the NSI couplings in three matrices according to 
whether the charged fermion in the external line is e, u, d 

(r-.e,u,d P p.e,u,d P p.e,u,d P 
t ee fc e/i t er 

p.*e,u,d P p.e,u,d P ^e,u,d P 
F *e,u,d P F *e,u,d P p e,u,d P 
fc er fc /iT fc TT 

Each set of three couplings e^g' d P enters in equation (11) in the entry v a p of the interaction 
potential matrix. Altogether there are 18 couplings with 36 parameters: each matrix of 
the three in eq.(18) contains 6 independent entries, each with a modulus and a phase. We 
analyse one coupling at a time, by taking all others zero. We first consider the cases of 
purely real and imaginary couplings, hence 4 possibilities for each phase 

0^ d = O,7r/2,7r,(3/2)vr. (19) 

Motivated by the arguments expound in the introduction we investigate the parameter range 
|e Q/J | e [5 x lCT 5 , 5 x 1CT 2 ]. We find that 

• Off diagonal entries e e ^ d P (a ^ (5) which contain 3x3x4=36 possibilities for moduli 
and phases do not induce any change in the LMA probability, nor any visible change 
in the rates, either if one or more at a time are inserted. 

• Diagonal entries e^ ,d p . 

(a) Real couplings e e ^ d p = ±\e e ^ d p \ (3x3x2=18 possibilities) do not change the 
LMA probability, hence the rates. 

(b) Imaginary couplings e^' d p = ±.i\e e ^ d p \ (3x3x2=18 possibilities) lead to proba- 
bilities which diverge from Plma for all |e aQ | > 5 x 10~ 5 . According to the probability 
shape that is obtained, we group these cases in the following way 
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Ta&Ze J - The NSI couplings that modify the LMA probability . 

All cases in the first column of table I along with D2, E2, F2 in the second column lead 
qualitatively to the same monotonically decreasing probability curve: a high probability 
(P > Plma) for low energy (E < 3MeV) and a low one (P < Plma) for intermediate and 
high energies. The curve becomes increasingly flat in this energy sector as e aa increases, 
which is also reflected in the flatness of the SuperKamiokande spectral rate. However for 
cases Al, Bl, CI and D2, E2, F2 the probability gets too high for low energies so that the 
Ga [4"5| |4*6] rate fails to be conveniently fitted. The 'best' results in the sense that they lead 
to the most flat spectral rate which approaches the SuperKamiokande one and to a correct 
fit for Ga are obtained alternatively from cases Dl, El, or Fl for the following values 

-z| £ f p | = -i 1.5 x 1(T 3 

-i\e u ^ p \ = -i 2.5 x 1(T 3 (20) 
+i\eif\ = +i 2.0 x 10- 3 . 

For larger values of the NSI couplings the probability moves further away from its LMA 
profile so that the Ga rate becomes too high and the 8 B one too low. In fig.l we plot four 
survival probabilities: at the lowest energy and from bottom to top the first curve is the LMA 
one, the next corresponds to all cases in (20) and leads to the best fit of the four, the next one 
to the case — i|eLf| = —i 3x 10 -3 and the top one to +i\el e p \ = +i 5 x 10 -3 . A comparison 
is shown in table II between the predictions and the quality of the fits obtained from LMA 
and the case — «|e^ p | = — i 1.5 x 10 -3 . In figs. 2 and 3 we show the SuperKamiokande 
spectrum for LMA (upper curves) and for the first case in (20) (lower curves) superimposed 
on the data points taken respectively from refs. [29J and [30]. The improvement obtained 
through the NSI coupling is clearly visible. In fig. 4 the two curves are superimposed on the 
SNO data points for electron scattering [47] . Here the data are also clearly consistent with 
a constant rate. 



Ga 


CI 


SK 


SNO NC 


SNO cc 


SNO E s 


Urates 


XsK sp 


XSNO 




LMA 64.9 


2.84 


2.40 


5.47 


1.79 


2.37 


0.67 


42.0 


48.6 


91.3 


-*Ku P \ 69.7 


2.74 


2.23 


5.47 


1.68 


2.26 


0.11 


40.3 


45.0 


85.4 



Table II - Comparison between the LMA predictions for solar event rates and the NSI ones 
with —i\s e ^ p \ = —i 1.5 x 10~ 3 . For details of the x 2 analysis see for instance [27^ . 
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The first set of cases in the second column of table I, namely A2, B2, C2, lead qualita- 
tively to the inverse behaviour with energy of the LMA probability. As \s aa \ increases from 
its lower bound, one gets P < Plma for low energies (E < 2 — 3M eV) and P > Plma for 
intermediate and high energies, so that the fits worsen with respect to the LMA ones. 

Finally all cases in the third column of table I, namely A3, B3, C3, D3, E3, F3 lead to 
probability curves which are totally unsuitable: they deviate drastically from both Plma 
and a from flat, suitable profile able to generate the SuperKamiokande spectrum. 

We have also checked that combinations of real and imaginary parts for all couplings do 
not change the previous results. This should be expected since, as mentioned earlier, purely 
real couplings do not change the LMA probability. The only consequence of introducing 
real parts in the NSI couplings comes in the spectral event rates through the quantities c/l,r 
in eqs.(14),(15), but the differences lie much beyond the experimental accuracy. 

3.2 {e%)cc±{s%)nc 

The analysis of the more general case of different CC and NC couplings affecting the v e e 
scattering diagrams can be made quite simple if one examines the coefficients of (s^)cc 
and {s^NC in eq.(ll). The first is unity whereas for equal CC and NC couplings it is 0.96 
and the second is now -0.04 as compared to the previous value 0.96 as well. Consequently 
one expects that the analysis for (£^)cc leads to approximately the same results as for 
equal couplings while the results are modified by a factor of 0.96/(-0.04) in the analysis for 
{ £ %)nc- Indeed the convenient modification in the LMA probability is obtained for 

-2(l^fl)cc = -2l.4xKT 3 (21) 

or alternatively 

+ *(l^ P |W = +*3.6xlO- 2 . (22) 

which lead to the same probability as the cases listed in eq.(20). The results for the other 
couplings involving u and d quarks are of course unchanged. As before we have considered 
one coupling at a time to be non zero. 

4 Conclusions 

We have investigated the prospects for improving the LMA predictions for solar neutrino 
event rates with NSI. At present there is no evidence of any new physics associated to a scale 
not too far above the electroweak scale, hence the great variety of theoretical models available 
for NSI. In our approach we assumed that NSI are extra contributions to the vertices v a vp 
and u a e, so the new couplings describe the deviation from the standard model. With this in 
mind we derived the neutrino interaction potential in solar matter which was added to the 
standard Hamiltonian, proceeding with the integration of the evolution equation through the 
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Runge-Kutta method. Neutral and charged current couplings are involved in interactions 
with electrons whereas only neutral couplings affect those with quarks. We considered the 
new interactions both at the propagation and at the detection level. The improvement we 
searched for the LMA predictions consisted in finding whether and how the modification 
induced by NSI can lead to a flat spectral event rate for SuperKamiokande and an event 
rate for the CI experiment within la of the data, while keeping the accuracy of all other 
predictions. 

We used the current notation for the NSI couplings e e a ^ d P where a, /3 are the neutrino 
labels and e, u, d denote the charged fermion involved in the process. We investigated the 
range \e a p\ e [5 x 1CT 5 , 5 x 1CT 2 ] and may summarize our main results as follows 

• Real couplings e e ^' d P = ±\e e ^ ,d P \ do not change the LMA probability when consid- 
ered either one at a time or altogether and thus they induce a small change in the 
neutrino electron scattering rate (< 1%) which is far beyond experimental visibility. 

• Off diagonal couplings e e ^' d P (a ^ (3) considered either one at a time or altogether 
do not change the LMA probability and thus the rates in a significant way. 

• Diagonal, imaginary couplings e^' d p = ±i\e^' d p \ are the only ones that lead to 
changes of all kinds in the probability and hence the rates. 

• The couplings that lead to the 'best' results in the sense of leading to a flat spectral 
SuperKamiokande rate and a good fit for all other rates follows 

-% 1.5 x 1(T 3 

-% 2.5 x 1(T 3 (23) 
+i 2.0 x 1(T 3 

within the assumption of equal charged (CC) and neutral current (NC) couplings. 

• If one assumes different (CC) and (NC) couplings we have 

-KK P \)cc = -zl.4xl0- 3 

+^(K P \)nc = +*3.6xl(T 2 . ^ 

In both cases (23) and (24) the results should be taken alternatively. 

Since real couplings do not change the LMA probability, the insertion of real parts in 
the couplings does not lead to any change in the above results. We also note that whereas 
an increase in sin9\^ leads to a mere decrease in the 8 B spectral rate which translates in 
a parallel shift of the curve [27], a convenient choice of NSI couplings [eqs.(23) or (24)] 
induces in contrast an increased flatness. The situation is similar to the magnetic moment 
conversion to sterile neutrinos [27] . 

Unless stated otherwise our analyses were done for one coupling at a time. There are 
of course infinite combinations of couplings so that a complete investigation is not possible 



— i\e 
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and would unlikely lead to new conclusions. The aim of this paper is to prove that NSI 
provide a true possibility to solve the remaining inconsistencies within the solar neutrino 
problem. 
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Figure 1: Survival probabilities as a function of neutrino energy in eV. At the lowest energy 
and from bottom to top, the first curve is the LMA one, the next is the one providing the best 
fit to the data with — i\e e ^\ = —i 1.5 x 1CT 3 [eq.(20)J, the next is for — i\s e ^\ = — i 3 x 1CT 3 
and the top one is for — f | = +i 5 x 10~ 3 with other non standard couplings vanishing 
in each case. The last two curves lead to an unacceptably high Ga rate prediction fjty \4^i - 
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Figure 2: Predictions for SuperKamiokande (units in Me V for electron energy). The upper 
curve is the LMA spectrum and the lower curve is the LMA spectrum with non standard 
interactions as in eqs.(23) or (24)- These are superimposed on the data published by the 
Collaboration in 2002 JEM/- 
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Figure 3: Same as fig .2 with the data published in 2008 [30]. 
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Figure 4: Predictions for SNO neutrino electron scattering superimposed on the data fj 7 ^ 
(units in MeV for electron kinetic energy). The upper curve is the LMA prediction and the 
lower curve is the LMA one with non standard interactions. Errors bars are larger than in 
SuperKamiokande so that the data are consistent with a flat spectrum. 
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